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^ , Abstract 

■ Central exclusive Higgs boson production, pp — > p(B H © p, at the LHC and s-channel 

resonant Higgs production in the photon-collider option of the ILC can provide a very im- 
portant contribution to the comprehensive study of the Higgs sector. Especially attractive 
is the bb Higgs decay mode, which for certain MSSM scenarios may become the discovery 
channel in exclusive Higgs production at the LHC and the Photon Collider (PC). Strongly 
suppressed and controllable backgrounds is an obvious requirement for the success of these 
exclusive measurements. One of the main sources of background comes from additional 
gluon radiation which leads to a three-jet bbg final state. We perform an explicit cal- 
culation of the subprocesses gg — > qqg and 77 — > qqg, where the incoming particles are 
required to be in a J z = state and the two gluons form a colour singlet, and investigate 
the salient properties of these potentially important background processes. 



1 Introduction 



The identification of the Higgs boson(s) is one of the main goals of the LHC Once the Higgs 
boson is discovered, it will be of primary importance to determine its spin and parity, and to 
measure precisely the mass, width and couplings. A comprehensive study of the whole Higgs 
sector, including precision mass and coupling measurements, spin and CP properties, will be the 
next stage. The conventional strategy to achieve this ambitious programme requires an intensive 
interplay between the LHC and the ILC (high-energy linear e + e~ collider) P . The ILC would 
enable a comprehensive phenomenological profile of the Higgs sector to be obtained, see for 
example j2J. In particular, a unique possibility to produce neutral Higgs bosons exclusively as 
s-channel resonances is offered by the 77 Compton Collider option of the ILC, see for example 
[HI IH EJ Ej and references therein. Whilst awaiting the arrival of the ILC, there has been 
growing interest in recent years in the possibility to complement the standard LHC physics 
menu by adding forward proton taggers to the CMS and ATLAS experiments (see for example 
|H1 El HOI HU [T2l HBl H31 US! Uni HZ! and references therein) . 




Figure 1: Schematic diagram for central exclusive Higgs production at the LHC, pp — > p+H+p. 



While experimentally challenging, this would provide an exceptionally clean environment to 
search for, and to identify the nature of, the new objects at the LHC. One of the key theoretical 
motivations behind these recent proposals is the study of so-called 'central exclusive' Higgs 
boson production, pp — > p®H®p. The © signs are used to denote the presence of large rapidity 
gaps; here we will simply describe such processes as 'central exclusive', with 'double-diffractive' 
production being implied. In these exclusive processes there is no hadronic activity between 
the outgoing protons and the decay products of the central (Higgs) system. The predictions 
for exclusive production are obtained by calculating the diagram of Fig. ^ using perturbative 
QCD [TBI EHIH]. m addition, we have to calculate and include the probability that the rapidity 
gaps are not populated by secondary hadrons from the underlying event [23 EE] • 

There are three major reasons why central exclusive production is so attractive for Higgs 
studies. First, if the outgoing protons remain intact and scatter through small angles then, to a 
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very good approximation, the primary active di-gluon system obeys a J z = 0, CP-even selection 
rule [221 122] • Here J z is the projection of the total angular momentum along the proton beam 
axis. This selection rule readily permits a clean determination of the quantum numbers of 
the observed Higgs resonance which will be dominantly produced in a scalar state. Secondly, 
because the process is exclusive, the energy loss of the outgoing protons is directly related to 
the mass of the central system, allowing a potentially excellent mass resolution, irrespective of 
the decay mode of the produced particle. 1 And, thirdly, a signal-to-background ratio of order 
1 (or even better) is achievable. As discussed in [23] , central exclusive production would enable 
a unique signature for the MSSM Higgs sector, in particular allowing the direct measurement 
of the Hbb Yukawa coupling. Moreover, in some MSSM scenarios this mechanism provides 
an opportunity for lineshape analysing ^3 H2j , and offers a way for direct observation of a 
CP- violating signal in the Higgs sector j2Hl C2] • 

The analysis in |TH1 Ell El was focused primarily on light SM and MSSM Higgs production, 
with the Higgs decaying to 2 b— jets. The potentially copious b— jet (QCD) background is 
controlled by a combination of the J z = selection rule [221 12"3] . which strongly suppresses 
leading-order bb production, colour and spin factors and the mass resolution from the forward 
proton detectors. It is the possibility to observe directly the dominant bb decay mode of the SM 
Higgs with Mjj ^ 140 GeV that first attracted attention to exclusive production at the LHC. 
It was subsequently realised that certain regions of MSSM parameter space can be especially 
'proton tagging friendly'. For example, at large tan/5 and M H < 250 GeV the situation 
becomes exceptionally favourable, with predicted Higgs signal-to-background ratios in excess 
of 20 ^7H2H]- In this particular case the tagged proton mode may well be the discovery channel. 
Though from an experimental perspective the bb channel is more challenging than the WW 
decay mode (see (2H1 EH [23 ), its unique advantages definitely merits a detailed analysis in 
realistic experimental conditions at the LHC. 

At the same time, the PC is especially best suited to the precise measurement of the 
T(H — > 77) width. Moreover, for certain regions of the MSSM parameter space, for example 
at the so-called 'LHC wedge', the PC has a discovery potential for the heavy pseudoscalar 
and scalar bosons, A and H, see for example 0. It is instructive to recall that in the case 
of 77 — > H — > bb production the potentially copious continuum b— jet background can be 
controlled by using polarised photon beams in the J z = initial-state (see for example |3Dl l3j). 
the same configuration of incoming particle polarisations that 'automatically' appears in the 
case of the pQCD box diagram of Fig. E for forward going protons at the LHC. The reason 
is that the Higgs signal is produced by photons (gluons) in a J z = state whereas the LO 
backgrounds are primarily initiated by the initial states with \ J Z \ =2, the J z = contribution 
being suppressed for large angles by a factor mf/s, see for example [211 130] . As discussed 
in [321 ESI f° r the 77 case, the physical origin of this suppression is related to the symmetry 
properties of the Born helicity amplitudes M^"** describing the binary background process 

7(Ai,/ci) + 7(^2,^2) -> q(Xq,p) + g(A ? ,p). (1) 

Here Aj labels the helicities of the incoming photons, and X q and Xq are the (doubled) helicities 

1 Recent studies suggest that the missing mass resolution a will be of order 1% for a 140 GeV Higgs, 
assuming both protons are detected at 420m from the interaction point |^] . 
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of the produced quark and antiquark. The fc's and p's denote the particle four-momenta, with 
s = (k 1 + k 2 ) 2 . 

Specifically, for a J z = initial state (Ai = A2) the Born quark helicity conserving (QHC) 
amplitude with Xq = —X q vanishes, 

M X x T Xq = 0, (2) 

see also [HI]. For the quark helicity non-conserving (QHNC) amplitude for large angle produc- 
tion we have 

M X x f" ~ o(^j M A A !!~ A % (3) 

where the amplitude on the right-hand-side displays the dominant helicity configuration of the 
background process. The above-mentioned m 2 suppression of the 77(</ 2 = 0) — > bb Born cross 
section is a consequence of Eqs. © and (J3J). The same is valid for the leading order amplitude 
of the gg pp — > bb subprocess, where the notation gg pp indicates that each active gluon in Fig. ^ 
comes from a colour- singlet t— channel (Pomeron) exchange and that the colour singlet di-gluon 
subprocess obeys the J z = 0, parity-even selection rule. 2 

The m\js suppression is especially critical in controlling the bb background. However, as 
was pointed out in [33], the suppression of the J z = background cross section is removed by 
the presence of an additional gluon in the final state. The radiative three-jet processes can 
then mimic the two-jet events in the quasi-collinear configurations when the gluon is radiated 
close to the b— quark directions 3 or (in the gg pp case) the extra gluon goes unobserved in the 
direction of a forward proton. First evaluations of the NLO QCD radiative backgrounds at the 
PC were performed in |33| I31] (for further development see (HlETj and references therein). This 
background contribution appears to strongly exceed the LO expectation and results in different 
shapes for various distributions. The background situation for the central exclusive H — > bb 
production at the LHC is much more complicated and requires a detailed combined study of 
various effects, see [TTH 114j. The analysis of some of these phenomena is still incomplete and 
require further detailed theoretical efforts, see, in particular, Sections 2.3 and 3 below. 

An important ingredient to this complex study which has not been completed so far is the 
availability 4 of the analytical expression for the matrix elements of the NLO process gg pp — > bbg 
which are needed to perform explicit calculations of the radiative background in the presence 
of realistic experimental cuts and selections. It is one of the main aims of this paper to derive 
the analytical expressions for the radiative cross sections, which then can be convoluted with 
existing Monte Carlo codes [3H1 E23 for the calculation of central exclusive processes. Note 
that, technically, the calculations of both gg pp and 77 induced colour singlet processes are 

2 It is worth noting that in the massless limit Eq. Q holds for any colour state of initial gluons. This is a 
consequence of the general property that the non-zero massless tree-level amplitudes should contain at least two 
positive or two negative helicity states, see for example |35| . It is an example of the more general Maximally 
Helicity Violating amplitude (MHV) rule, reviewed for example in |36j . 

3 In the PC case there also a sizeable radiative background coming from ccg production. 

4 Although a number of automated packages are available for tree-level scattering amplitudes for arbitrary 
final states, it is very difficult to extract from these the projection onto a specific spin (e.g. J z — 0) and colour 
(e.g. colour singlet) initial state. 
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quite similar, the latter providing just a subset of the diagrams for the former. Therefore, it is 
convenient to discuss the two radiative processes simultaneously, illuminating their similarities 
and differences. These are discussed in Sections 3 and 4. 

Our phenomenological discussion will be focused on central exclusive bb production at the 
LHC. However it is also worth noting that the CDF Collaboration at the Fermilab Tevatron 
is currently performing an experimental study of exclusive diffractive bb events |HJ| 5 and the 
results of this paper could prove useful for the analysis of these measurements. 

2 On the backgrounds to the p+ (H — > bb) +p signal 

2.1 Classification of the backgrounds to exclusive Higgs production 

From the theoretical point of view, it is convenient to consider separately the QHC and QHNC 
amplitudes. These amplitudes do not interfere, and their contributions can be treated inde- 
pendently. This could be especially convenient at the stage when the parton shower algorithm 
is applied, since double counting can be avoided. 

There are two main sources of gg pp — > bb background process: 

(i) the LO 0{a%) QHNC amplitude squared, 

(ii) the NNLO O{og) QHC contribution which comes from the one-loop box diagrams. 

As already mentioned, there is also the possibility of a NLO 0(a%) gg pp — > bbg background 
contribution, because large-angle, hard-gluon radiation does not obey the selection rules. Of 
course, the extra gluon may be observed experimentally in the central detector, and such 
background events eliminated. However, there are important exceptions which we discuss 
below. 

In the case of the NLO gg pp —>■ bbg process the dominant contribution comes from the QHC 
amplitude, since the QHNC piece is additionally mass-suppressed. Here we consider two types 
of radiative background process that can mimic the H — > bb central exclusive signal. 

(a) The extra gluon may go unobserved in the direction of one of the forward protons. This 
background may be reduced by requiring the approximate equality M missing = M bb . But 
the degree of this reduction will depend on the mass resolution in the proton detector and 
jet energy resolution in the central detector. Since the mass (jet energy) resolution AM bb 
in the central detector is expected to be much worse than the missing mass resolution, 
AM missing <C AMjj, the background will be limited in practice by the AM bb value. 

(b) The remaining danger is large-angle hard gluon emission which is collinear with either 
the b or b jet, and, therefore, unobservable. As discussed in (TUJ EHj, f° r Jz — this is 
suppressed for soft gluon radiation. Although there is a certain suppression of collinear 
radiation as well, this issue requires further detailed analysis, see below. 

According to the study in [TO] , if the cone angle needed to separate the g jet from the 
b (or b) jet is AR ~ 0.5 then the expected background from unresolved three jet events 
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leads to B/S — 0.2. The calculations presented in Section 3 below will allow a more 
precise evaluation of this ratio. 

It is worth noting that a detailed experimental study of the three-jet bbg final state 
could be useful for background calibration purposes. This is of particular interest for 
the kinematic configurations which are enhanced, for example when an energetic gluon 
recoils against a quasi-collinear bb pair. Recall also that the CDF collaboration is currently 
measuring exclusive bb production at the Tevatron |4Uj . 

Note that in this paper we do not discuss the effects coming from collisions of two soft 
Pomerons, neither do we address a possible contribution from central inelastic production, see 
[§| . The reduction of such backgrounds is controlled by imposing the missing mass equality, see 
|10j . Note also that gluon radiation off the screening gluon (labelled C Q' in Fig.[IJ is numerically 
small jH]. 

2.2 Properties of the leading-order gg pp — > bb background process 

As we have discussed, an important advantage of the p + (H — > bb) + p signal is that there 
exists a J z = selection rule, which requires the LO gg pp — > bb background to vanish in the 
limit of massless quarks and forward outgoing protons. However, in practice, LO background 
contributions remain, see The prolific gg pp — > gg subprocess can mimic bb production 

when the outgoing gluons are misidentified as b and b jets. Assuming the expected 1% proba- 
bility of misidentification, and applying a 60° < 9 < 120° jet cut, gives a background-to-signal 
ratio B/S ~ 0.2 [Ej. (Here and in what follows, we assume for reference that the mass 
window over which we collect the signal is AM ~ 3a = 3 GeV). 5 

Secondly, there is an admixture of \J Z \ = 2 production, arising from non- forward going 
protons, which gives B/S ~ 0.05, see Section 2. 5. 6 

Thirdly, in reality the quarks have non- zero mass and there is a contribution to the J z = 
cross section of order ml/Ej,, where Et is the transverse energy of the b and b jets. In [TU1 E] 
the contribution from this source was estimated as B/S ~ 0.2. However the higher-order 
QCD effects may strongly affect this result. First, there is a reduction coming from the self- 
energy insertion into the fe-quark propagator, that is from the running of the b— quark mass 
from mb(mb) to its value fn^Mji) < rf%b(mb) at the Higgs scale. Here mb(n) is the running 
b— quark mass in the MS scheme |42j . It is known that in the H — > bb decay width these 
single logarithmic (SL) (a^ln^ 2 -) effects diminish the corresponding Born result by a factor 
of approximately two |4"3*| . Although this still requires a more formal proof, we strongly believe 
that the same (factor of two) reduction applies in the case of mass-suppressed J z = binary 
reactions gg pp — > bb and 77 — > bb at large angles. 

There is another (potentially important) source of uncertainties in the evaluation of the 
rate of exclusive bb production at J z = 0. This is related to the so-called non-Sudakov form 

5 Such a background is practically negligible in the case of the PC since it must be mediated by the higher- 
order 'box' diagrams. 

6 Analogous to this in the PC case is the contribution from the initial photon state with | J z \ = 2, which may 
constitute a non-negligible source of background. 
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factor in the cross section F q which arises from virtual diagrams with gluon exchange between 
the final quarks (or initial gluons), see m the 77 case the double logarithmic (DL) 

approximation to F q has the form 

F q (L m ) = £c n (^) n (4) 

V 7T J 

with 

L m = In , (5) 

Co = 1 and c\ = —8 03] so that the second (negative) term in is anomalously large and 
dominates over the Born term for Mh ~ 100 GeV. This dominance undermines the results of 
any analysis based on the one-loop approximation. The physical origin of this non-Sudakov 
form factor was elucidated in |33] where its explicit calculation in the two-loop approximation 
was performed. It was also shown that for reliable calculations of the DL effects the two-loop 
calculation should be sufficient. This was subsequently confirmed by a more comprehensive 
all-orders study [33] ■ As is well known, there are other DL effects (the so-called Sudakov 
logarithms j3E]) that arise from virtual soft gluon exchanges. Their contribution depends on 
the particular kinematics in the final state. As discussed in [33], in the case of quasi-two-jet 
configurations Sudakov and non-Sudakov effects can be with good accuracy factorised, because 
they correspond to very different virtualities of the internal quark and gluon lines. For the 
final state radiation, the Sudakov effects (both for the signal and for the background) can be 
implemented in parton shower Monte Carlo models in a standard way. For the gg pp initial 
state, the Sudakov factors are explicitly incorporated in the unintegrated gluon densities, see 
[HI HE]- Currently, for Higgs production at the PC the DL factors are accounted for by the 
simulation programme used for generating background events, see [37J. Unfortunately, from a 
phenomenological perspective, it seems to be potentially dangerous to rely on the DL results, 
since experience shows that formally subleading SL corrections may be numerically important. 
We plan to address this issue in future. 

Nevertheless, as a first step in understanding the situation in the pp case it is instructive to 
evaluate the size of the DL effects for the gg pp — > bbg reaction. Recall that in the photon-photon 
case the two-loop expression for F q takes the form [33] 

F q = (l-3^ + ^ (l + ^), (6) 

with 7 

T = ^ C F L 2 m . (7) 

7T 

The corresponding one-loop result for the gg pp initiated process is 

F g = 1 - (2CV + 4N C )F/C F ~ (1 - (C F + 2N C )F/C F ) 2 (8) 

7 Note that an additional problem is that while in the Born cross section it is natural to evaluate as at the 
hard scale Mh, for T we have no reason to adopt this prescription. The existing PC generators do not take 
into account possible differences in the scale of as for the different quantities. 
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Unfortunately, due to the large colour coefficients the one-loop DL contribution becomes 
larger than the Born term, and the final result will be strongly dependent on the NNLO 
contribution as well as on the scale \i at which the QCD coupling as is evaluated and on the 
running b— quark mass. It seems plausible to choose the scale n ~ M#/2 (since we are interested 
in the bb background at s = = Mjj) for the Born amplitude, but a lower scale /i ~ y/Mnrrib 
for the factor T in Eq. (J7|) which originates in the region where the quark propagators are close 
to the mass-shell. The NLL can be effectively incorporated in Eqs. (I6I8|) by introducing a scale 
factor c in the argument of the logarithm, that is by replacing the ratio M#/m& by cMn/mb. 
It follows from the comparison with the complete one-loop calculation jH] for the process 
ll{Jz = 0) — > bb, that the scale factor c ~ 0.5. This looks quite reasonable if we account for 
the kinematical configuration. 

Without the complete result for the higher-order radiative corrections corresponding to the 
gg pp — > bb amplitude, it is impossible to make a firm prediction. To gain an insight into the 
size of the possible effects, we make the assumption that the same scale factor c = 1/2 is 
valid in this case as well. Then choosing as = a^M^m;,) ~ 0.15, the value of the correction 
(1 + 2N c /Cf)J 1 ' ~ 2.5 exceeds the Born term. In other words, the whole amplitude changes 
sign and the background cross section becomes a few times larger than the Born expectation. 

Accounting for the running 6-quark mass, the expected non-Sudakov correction factor can 
be approximated by 



■ _ MMnnib) a s (M H m b ) 2 / c - M H \ 

m b {M 2 H ) cJ 7T \m b (M H m b )J 



(9) 



For the 77 case we have a similar expression 



_ m b (M H m b ) a s (M H m b ) 2 / c ■ M H \ 
m h (M 2 H ) F it {m b (M H m b )) 



(10) 



If we take these formulae literally we would conclude that in the gg case the quasi-two jet cross 
section is a factor of 2 larger than the 'naive' (but frequently used) Born prediction, calculated 
with as(M H ) and the b— quark pole mass. Similarly, for the 77 process the result is about 5 
times lower than such a naive Born estimate. 

While for the 77 case the estimated effect seems to be reasonably justified, for the gluon- 
initiated process it can serve only as a rough illustration of the possible size of the effect. The 
actual results will depend crucially on the value of the scale factor c (i.e. the NLL contributions) 
and on the specific choice of the arguments of the running coupling and 6-quark mass. The 
main purpose of the exercise was to demonstrate that, currently (or, at least, before a complete 
one-loop result becomes available), the gg pp — > bb cross section can be estimated to no better 
than an order of magnitude accuracy. 



2.3 Quasi-two-jet-like radiative background events 

As was first found in |H] for the 77 case, there is an additional NNLO contribution which is 
not mass-suppressed and is potentially important especially for large energies. It comes from 
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the QHC box diagrams. This piece cannot be evaluated in terms of the tree-level amplitudes 
using the cutting rules in 4 dimensions. A systematic method for calculating such amplitudes 
in the massless limit is based on generalised unitarity in .D- dimensions, see for instance |47j . 
An explicit calculation in [44] using dimensional regularization of the NNLO process 77 — > qq 
at large angle 6 gives 

^ (ll^qq,J z = 0) = ^C 2 ^ cos 2 6(1 -cos 2 6) (11) 

with 

da Bm _ lWg* (3(1 -/?*) 

^(77 - 11, J, = 0) = — {1 _^ CQs2er (12) 



where (3 = yl — Am 2 / s, and m q and Q g are the mass and electric charge of the quark respec- 
tively. 8 Note that as is easily seen from Eq. (|12j) . the NNLO elastic cross section vanishes at 
6 = 7r/2. This is a consequence of the rotational invariance about the quark direction at 180° 
and the identity of the photons (see Ref. [33]) and remains valid in the absence of radiation 
at all orders in as- This is also true for the gg pp — > qq process. The ratio in ((TJ)) reaches its 
maximum at 6 = 7r/4 where 

foNNLO 2 

(77 - qq, J z = 0) = ^~ 2 — 2 . (13) 

da Born 72tH ml 

Accounting for the running b— quark mass and the NLLO effects discussed in the previous 
Section, we conclude that even at M# ~ 140 GeV these NNLO contribution to the cross 
section could not exceed 0.1 of the modified Born term. The NNLO elastic cc contribution is 
16 times larger, and at Mh — 130 GeV becomes comparable with the modified bb exclusive 
term. However with a reasonably good experimental c— quark rejection this background can be 
strongly reduced without seriously degrading the H — ► bb signal. 

Using the existing results for the one-loop amplitudes of the gg —> bb process in the massless 
limit (see for example ^Hl[49]) we can write down the corresponding NNLO expression for the 
ratio of the gg pp — > bb subprocesses as 



da NNLO . PP ^ (C F -N C ) 2 da NNLO , 

{gg -> qq, h = o) = K n2 1 ■ -, (77 -> qq, J* = 0). (14) 

«C Born <~<F Ct<J Born 

Note that the appearance of the (Cp — N c ) 2 factor in Eq. (|T4*j) is not accidental. It is a con- 
sequence of supersymmetry requiring the vanishing of such helicity amplitudes in a supersym- 
metric theory, which happens if we put the fermions in the adjoint representation (gluinos). 9 
This is in marked contrast with the combination (Cp+2N C ) 2 that appears in Eq. (|5jl. where the 
result is of a purely classical nature and supersymmery arguments cannot be applied. Note that 
in the massive quark case, even if we were to consider altering its colour representation from 



8 In |34| there are some confusing statements regarding the properties of the 77 — > bb amplitude in the 
complex plane and the interpretation of the one-loop amplitude result of 01] . But these do not affect the actual 
formulae. 

9 We are grateful to Lance Dixon for an illuminating discussion of the properties of helicity amplitudes in a 
supersymmetric theory. 
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the fundamental to the adjoint representation, we cannot put it into the same supersymmetric 
multiplet with the massless gluon [50]. 



2.4 Admixture of | J z \ = 2 production, arising from non-forward go- 
ing protons 

In the exact forward direction, the J z = selection rule is just a consequence of the s-channel 
helicity conservation for the forward protons. For the non-zero transverse momenta of the 
outgoing protons (pi t t,P2,t) some admixture of the \J Z \ = 2 component appears. Its value is 
controlled by the orbital momentum transfered through the Pomeron (i.e. through the colour 
singlet t-channel two gluon exchange) that is by the product {j> t) i 4 r t), where i = 1,2, and the 
effective Pomeron size r t ~ 1/Qt is driven by the inverse transverse momentum Q t in the gluon 
loop. Thus, this admixture of the \J Z \ = 2 states can be evaluated as j22] 2pi,tP2,t/Qt in the 
amplitude or 

— o?— < 15 > 

for the cross section. An additional factor 2 arises from the azimuthal angular averaging 
{Q 3 Q k ^5flQV2). 

Technically the polarisation structure is as follows. In the equivalent gluon approximation 
the polarisation vector of the active gluon is proportional to its transverse momentum e Mj j ~ 
(pt,i ~ (~^YQt)fi/xi (i = 1,2). 10 Thus the product of two polarisation vectors can be written 

as 

(pm - Qt)»(pt,2 + Qt) v 

eu,ie„2 — , 

XiX 2 

e^ie Ut2 oc p t ,i^Pt,2,u - Qt,nQt,v + (Qt 7t iPt,2,u - Pt,utQt,u)- (16) 

After the (Q t ) angular integration the last term in (|16j) vanishes while the second term gives 
—8^Q\/2. In terms of helicity amplitudes, 5j?J corresponds to a pure J z = state. On 
the other hand, the first term in after the averaging over p t ^ directions, generates the 

J z = and \J Z \ = 2 states with equal probabilities. Unfortunately, the amplitude with an 
extra 1/Qf factor becomes less convergent at small Q\ . The dominant contribution comes from 
the region of relatively low Q t ~ GeV (and even lower for the Tevatron energies). Therefore, 
we cannot guarantee the precision of the numerical evaluation. Using the MRST99 partons 
[5T] we expect the \J Z \ = 2 admixture for central exclusive production of the state with the 
mass M ~ 120 — 160 GeV to be about 5% at Tevatron energies and ~ 1.5% at the LHC 
(y/s = 14 TeV) 123] . 

There is good news however. It is worth noting that the \J Z \ = 2 contribution to the 
gg pp — > qq background is additionally suppressed numerically (by a factor of, at least, ~ 0.2). 11 . 
In order to gain an insight into the origin of this additional suppression, we note that the cross 
section vanishes at 9 = tt/2 in the qq rest frame (neglecting the proton transverse momenta 

10 Alternatively, the same result can be obtained in the LO using the planar gauge (n^ • A^) = with the 
gauge 4-vector n M parallel to the 4-momentum of the centrally produced system M. 

11 This suppression was not accounted for in [2B|23|- This provides added value to the improvement of the 
signal-to-background situation in the 66-case. 
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in comparison with the transverse energy of the quarks). This follows from the identity of the 
incoming gluons (protons) and invariance with respect to the 180° rotation about the quark 
direction. This in turn, causes the cross section to be proportional to cos 2 9. This phenomenon 
is also seen in the vanishing of the NNLO non-radiative amplitude at 9 = n/2 considered in 
Section 2.3, and in the soft radiation off the screening gluon, considered in jH] (see also the 
discussion in fHBj). 

2.5 NLO radiation close to the beam directions 

The NLO subprocess gg pp — > bbg can also avoid the J z = selection rule. Extra gluon radiation 
in the beam direction goes into the beam pipe, and cannot be observed directly. Therefore, 
experimentally, the event may look like central exclusive production. There are two main 
consequences of this extra gluon radiation: a) the system M 8 which is centrally produced via 
gg — > Ms fusion is now in the colour octet state (which we label by the symbol '8'), and (b) 
the J z = selection rule, which suppresses the bb LO QCD production in the genuine central 
exclusive event, becomes redundant. 

Let us discuss this point in more detail. The emission of a low qt extra gluon is strongly 
suppressed due to the distructive interference between the amplitudes where the gluon q^ is 
emitted from the right (active) or from the left (screening) t-channel gluons in Fig.0 Therefore, 
it is sufficient to consider the case when Mh ^ qt ^ Qt, and the polarisation of the active gluon 
(participating in the gg — > M 8 fusion process) is directed along the new vector q t . Assuming 
that the extra gluon with momentum g M is emitted from the lower active gluon (i = 2), the 
polarisation structure becomes 

e^ie Vt2 oc (jh,i,h - Qt,n)(qt, v + Qt,u) — Pt,i,^qt,u - Qt,nQt,u - 6$Q 2 /2 , (17) 

where the last term corresponds to the J z = state and the first term corresponds to the hard 
subprocess, which after the q t and p t ,i averaging looks like the usual fusion of two unpolarised 
gluons (with equal probabilities for the J z = and \J Z \ = 2 initial states). The contribution 
generated by the second term is more complicated. At first sight it should vanish after the 
integration over the azimuthal Qt angle, however due to the factor (Qt + qt) 2 in the denominator 
of the amplitude some component of the momentum Q t)fi in the g t|jU direction still survives. In 
the limit of q t ^> Q t this leads to a contribution ~ Qt,nQt,uQt / Qt t° the right-hand side of 
Eq. ([T7|). which again contains the J z = and | J z \ = 2 states with equal probabilities. 

Thus, the ratio of the | J z \ = 2 to J z = contributions to the cross section for this process 
may be evaluated as 

*(l^l=2) <Q 2 t > 2 + <p 2 Xq 2 t > 

a(J z = 0) " 2<Q 2 >2 + <pf ><??>■ 1 ) 

We might expect that the probability to emit such an extra gluon would contain a double 
logarithm, but this does not happen. First, as will be discussed in Section for the massless 
b— quark case the soft gluon emission gg pp — > bbg is suppressed by a factor (Eg/M b i) A and 
therefore gives no logarithm. On the other hand, the collinear logarithm is limited by the 
angular acceptance of the detector. Any gluon with a sufficiently large q t will be observed 
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in the Central Detector. Moreover, up to pseudorapidities \r]\ ~ 6 — 7 the extra gluon jet 
will be observed experimentally in a Forward Detector. Such events with a third jet will be 
easily distinguished from the H — > bb decay. Next, if the energy of the gluon q exceeds the 
mass resolution then there will be no matching between the missing mass calculated from the 
momenta of the outgoing forward protons and the mass M$ measured in the Central Detector. 
Assuming the mass resolution AM 6 g ~ 20 GeV, we require that for central bb production the 
energy of the third (forward) gluon jet must be less than 40 GeV, and to get \rj\ > 6 such a jet 
must have a very small transverse momentun, q t < 2E g exp(—6) = 0.2 GeV. The production of 
such a low (?t(<C Qt) gluon is strongly suppressed by the interference between the emissions of 
the active gluon q t and the screening gluon Q (see Fig. HJ. This contribution becomes smaller 
than the admixture of the \J Z \ =2 states. 



3 The <7<7,77 — > qqg \J Z \ = colour singlet hard process 



In this section we will present results for the matrix elements squared for the colour singlet 
hard scattering processes gg — > qqg and 77 — > qqg. Since we will be using these results in 
situations where the momentum transferred in the hard scattering is much larger than the 
b— quark mass, we will set m q = 0. We will compare our results with the corresponding full 
spin- and colour-summed amplitudes in order to exhibit the different limiting behaviours. 

In fact the spin- and colour-summed matrix element squared for the 2 — > 3 process g(p\) + 
9(3*2) ~^ 9{Ps)q.{P4)<i(Pb) process has been known for a long time [52! and has a relatively simple 
analytic form: 



E\m\ 



9\ 



AN 2 (N 2 - 1) 



' aibi(a\ + b\) + a 2 b 2 (a\ + bj) + a 3 6 3 (a^ + b\ 
aia 2 a 3 &iM3 



- + N \ - - 



d 



12 



d 



23 



d 



13 



| 2iV 4 / a 3 ^3(Qi^2 + Q2&1) Qi&i(Q2&3 + ^3^2) 0^2(^1 + aib 3 ) 



<^23 ^13 



di? d 



12 "13 



dio d 



12 "23 



(19) 



where a { = p { ■ p 4 , b { = p { ■ p 5} d {j = p { ■ pj = 1, ...3) and s = 2p 4 ■ p b . An averaging 
over initial spins (1/4) and colours (l/(iV 2 — l) 2 = 1/64) has been performed. In fact the 
above spin-summed amplitude squared comprises 12 distinct non-zero helicity combinations, 4 
of which 

(++; - - +), (++; - + -), (— ; + - +), (— ; + + -) (20) 

(in an obvious notation) correspond to a J z = initial state, while the remaining 8 



-+; + -+), (-+; 
(+-; + -+), (+-; 



+), (-+; + + -), (-+;- + -), 
-+), (+-; + + -), (+-;- + -) 



(21) 



correspond to a | J z \ = 2 initial state. Note that in all cases X q = —Xq, corresponding to helicity 
conservation along the fermion line. The other important point to note is that all the above 
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combinations are MHV amplitudes [35], m the sense that the sum of the (five) helicities is 
always ±1. At this order there are no NMHV amplitudes or higher. MHV gggqq scattering 
amplitudes have a very simple analytic form, see the Appendix, and so with appropriate colour 
weightings and momentum permutations, the colour singlet J z = 0, ±2 matrix elements squared 
can be easily constructed. The result is 

2 4 

Y,\M\ 2 (J Z = 0; colour singlet) = - £ \z(l, 2, 3, h) + z(2, 1, 3, h) + z(3, 2, 1, h) 

9 h=l 

+z(3, 1, 2, h) - ~ (*(1, 3, 2, h) + z(2, 3, 1, h)) f (22) 

where spin and colour averaging factors are included. The z factors are given in the Appendix. 
Expressions for the other spin and colour combinations can also be written down in terms of 
the z(i,j, k, h) factors. However these are more lengthy and so will not be presented here. 

Compact analytic expressions exist for the corresponding 77 — > gqq spin summed and J z = 
amplitudes squared. In the notation of Eq. (fTTJj) . 

1.^12/ ■ in 2 4 s f a>ih(al + b\) + a 2 h 2 [a\ + &|) + a 3 b 3 (aj + bj)\ 
£|AC(SP- -nmed) = 8 ff ,V - { a^al^l ~ J • 

EW, = 0) = B^£(^±|). (23) 

4 Numerical results and discussion 

In this section we consider some of the properties of the amplitudes presented in the previous 
section and, in particular, focus on the differences between the spin-summed (unpolarised) and 
J z = cases. The differences are most dramatic for the kinematic configuration in which the 
final-state gluon is soft. It follows from the E g —>■ limit of Eq. (J23|) that the matrix element 
squared for the J z = case is proportional to Eg, while in the unpolarised case it exhibits the 
standard 1 /Eg behaviour. This is because the first two terms in the bracket on the right hand 
side in the unpolarised case, which are responsible for the leading infrared behaviour in the 
soft-gluon limit, are absent in the J z = case. The net difference is four powers of a 3 or 6 3 , 
equivalently Eg. Numerical calculation shows that exactly the same behaviour is found in the 
gg scattering cases. 

In terms of the cross sections in the soft gluon limit, 

dE g E " (24) 

while in the unpolarised case we arrive at the the standard infrared behaviour 

^""unpol 



OlEg Eg 



(25) 



Such behaviour is rooted in the Low-Burnett-Kroll (LBK) [SH] theorem (see also JH3EHE3)- 
According to the LBK theorem, for radiation of a soft gluon with energy fraction x g <C 1, 
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the radiative matrix element M ra d may be expanded in powers of the scaled gluon energy 
x 9 = Eg/Eb 

1 oo 

M rad ~ — C n x n g , (26) 

X 9 n=0 

where the first two terms, with coefficients Cq and C\ (which correspond to long-distance 
radiation), can be written in terms of the non-radiative matrix element Mb- The application of 
these classical results is especially transparent when the cross sections are integrated over the 
azimuthal angles. Then the non-radiative process depends only on simple variables, such as the 
centre-of-mass energy. When Mb = 0, the expansion starts from the non-universal C2X 2 term, 
which corresponds to non-classical (short- distance) effects, not related to Mb- This is exactly 
the case for the J z = Born amplitudes which vanish for massless quarks for both 77 and gg pp 
processes. 12 On the other hand, in the soft limit the unpolarised result is dominated by the 
non- vanishing non-radiative amplitudes, either ( — h; H — ), ( — h; — h) or (H — ; H — ), (-| — ; — h), 
and in this case the matrix element squared, |M rac j| 2 , is proportional to 1/-Ey Recall, however, 
that the mass-suppressed contributions to the J z = amplitudes will induce the normal infrared 
behaviour, see for example J33], i>e. 

MJ Z = 0) _ (27) 

(lEg ~ S Eg ^ } 

The above behaviour is illustrated in Fig. |2] where we plot 13 the matrix elements squared 
as a function of the scaled gluon energy x g = E g /E heam . For simplicity, we work in the centre- 
of-mass frame with all final state particles in the transverse plane and equal energies for the 
(massless) quark and antiquark. In this configuration, x g — > is the soft gluon limit, x g = 2/3 
is the (transverse) 'Mercedes' configuration, and x g — > 1 corresponds to a transverse gluon 
balanced by a collinear qq pair. 

The expected LBK behaviour discussed above [x 2 g versus x~ 2 ) is clearly seen in the x g — ► 
limit. Note also that the gg (but not the 77) amplitudes become singular as x g approaches the 
kinematic limit at x g — 1. This is the collinear singularity caused by the production of two 
final-state back-to-back gluons, one of which splits into a quark-antiquark pair, i.e. g — > qq. 
For these kinematics, this is manifest as a (1 — x g )~ l singularity in the matrix elements squared. 
There is no analogue in the 77 case, and indeed here the amplitudes vanish in the x g — > 1 limit, 
since at x g — > 1 the quark and antiquark go in exactly the same direction and their electric 
charges screen each other. Thus the coupling of the photon to the qq pair vanishes. 

More generally, we can study the behaviour of the matrix elements squared when all three 
particles lie in the transverse plane by using the Dalitz-plot variables x q and x g , the scaled 
energies of the final state quark and antiquark respectively, with x q + x g + x g = 2. For the 
J z = 77 — >■ qqg case, it is straightforward to show that 

E \MUJ Z = 0) = 2 *£* [(1 - x q f + (1 - x q f] X " + Jy\ (28) 

12 Actually, the J z — non-radiative matrix element vanishes for any colour state of two gluons since it 
corresponds to a Maximally Helicity Violating (MHV) situation, see for example 
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For display purposes, the QCD and QED couplings g s and e are set to 1 in this plot. 
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In the limit where Xg (or x q ) — > 1, an antiquark balances a collinear quark-gluon pair, and the 
matrix element squared becomes 

J2\M\^(J z = 0)^^^(l-xr/x, (29) 

with x = 1 — x g the fractional quark momentum of the quark-gluon pair. Note also that 
the matrix element squared vanishes when the quark and antiquark are emitted in the same 
direction (x q + x q = 1), as anticipated above. 

For the gg scattering case, there is no simple analogue of (|2*Hj) for the whole Dalitz plot. 
However, the behaviour along the boundaries can be extracted and studied. We find 

J2\M\ 2 gg (j z = o) = 9 -f {(i-x q y + (i-x q y] F{ *f^ (so) 

where F(x, y) is a non-singular function with constant values along the three sides of the Dalitz 
plot: 

64 

F = — , when x = 1 or y = 1 
9 y 

F = 36, when x + y — 1. (31) 

Note the jump of the function F at x = l,y — > and y = l,x — > 0. These are the points 
where the soft antiquark (or quark) changes its direction, leading to a discontinuity in the 
classical coloured current. 

To summarise, in the gg J z = colour-singlet case, the only final-state 14 singularity is 
when the gluon is emitted opposite a collinear quark- antiquark pair. In practical terms, this 
corresponds to the case when both 6-quarks are contained within the same jet. Strictly speaking, 
this does not constitute a background to Higgs production, since the latter gives rise to two 
distinct 6-jets. 

In contrast, the full spin-, colour-summed gg amplitude has additional singularities when 
x q — 1 or Xq — 1. In fact the full singularity structure is exhibited in the empirical form 

£ |AC(sP.n. -lour summed) = f _ , (32) 

where G is non-singular throughout the Dalitz plot. Note that when /2 — the 

kinematics of Fig. |21 — this reduces to 

J2 \M\ 2 Js P m, colour summed) = G{1 - xj2,l -x g /2) ^ ^ 

s x g\*- x g) 

What does this mean for jet cross sections? Recall that in e + e~ — > qqg the 'three-jet 
cross section' is defined by integrating the matrix element squared over a region away from the 



14 We note that we are considering here central production of all three final-state particles. There are of 
course additional singularities when the final-state gluon is emitted collinear with the incoming gluons, but 
these initial-state singularities are well understood. 
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Figure 2: Dependence of the gg,77 — > bbg matrix element squared on the final-state gluon 
fractional energy, when all three particles are produced in the transverse plane and the b and 
b have equal energy. Note that the b— quark mass is set to zero. 
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boundaries of the Dalitz plot, the exact region depending on the jet algorithm definition. For 
example, the y cnt (JADE) algorithm defines the three-jet region by 

1 - X q , 1 - Xq, X q + Xq + 1 > y cn f (34) 

Given the absence of collinear singularities in the J z = 0, colour- singlet case when the gluon is 
emitted parallel to the quark or antiquark, we may therefore expect relatively more radiative 
three-jet events than in the \J Z \ = 2 or spin-summed cases. Note, however, that the three 
jet cross section is still (logarithmically) singular in the y cut — > limit corresponding to the 
configuration x q + Xq — > 1 in which the final-state b and b are collinear, i.e. g+g — > g+g*{—* bb). 
However, such radiatively generated quasi-two-jet events can be suppressed by requiring two 
distinct, spatially-separated b— tagged jets. Indeed applying a minimum-angle cut between the 
b and b jets leads to a finite (J z = 0, colour-singlet) three-jet cross section. 

Note that in order to compare the matrix elements for the 77 — > bbg and gg — * bbg reactions 
in Fig. 2 we set the couplings g s = e = 1. Strictly speaking, at the leading order (tree level) 
at which we work, it is not known at what scale the couplings must be evaluated. However in 
actual cross section calculations it would appear natural to take two of the vertices at the large 
scale n\ 2 ~ s/2 and the third (outgoing gluon emission) vertex at a lower scale /i§ ~ /c^, i.e. 
\M\ 2 oc a%(s/2)a s (kj_), where k± is the transverse momentum of the final gluon with respect 
to the nearest (quark, antiquark or beam) jet direction. 

In summary, provided that two distinct b— jets are required, then the dominant background 
arising from gg — > bbg production in the J z = 0, colour-singlet case corresponds to three-jet 
production. The same is true for 77 production. Although we have concentrated our analysis 
on production in the transverse plane, this conclusion is still valid for central production, i.e. 
where the final state jets are restricted to a central region in rapidity. 



5 Summary 

In the previous sections we have shown(see also [^2]) that, as a consequence of the Low-Burnett- 
Kroll theorem when neglecting the quark mass, the differential distribution over the gluon 
energy in the 77, gg pp — > qqg cross section is 

dE g ~ E " (35) 

in marked contrast to the Higgs or unpolarised case, where the cross sections exhibit the 
standard infrared behaviour 

C^unpol 1 



Eg 



(36) 



As a result, the relative probability of the Mercedes-like configuration in the final qqg state for 
the J z = background processes becomes unusually large. We have derived explicit analytic 
expressions for the 77, gg pp — > qqg amplitudes using MHV techniques, and calcualted some 
simple energy distributions to illustrate their generic behaviour. These amplitudes can easily 
be incorporated into more sophisticated Monte Carlo programmes to investigate background 
event rates in the presence of realistic experimental cuts. 
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Finally, the approach of Ref. [SH] enables us to evaluate the difference between the charged 
multiplicities of the signal N$ and Mercedes-like background events N^q c containing 6-quarks 
for the J z = initial state for both processes 77 and gg pp induced processes. As is shown in 

m, 

AN = K-(M H ) - N S (M H ) = = Ngg + l -N gg (^) - N m (M H ). (37) 

For example, for a 100 GeV Higgs boson, 2E* = ^M- ~ 58 GeV, which corresponds to the 
energies of the existing measurements by TOPAZ and VENUS, see [56j. Substituting into 
Eq. (}3*Tj) the corresponding experimental results for N q q and the fits to the gg multiplicity from 
|56j . we arrive at the multiplicity difference between the Mercedes-like background events and 
the bb signal, 

AN =6.8 ±1.5. (38) 

A similar result (AN ~ 8.0) appears if we use the existing (udscb) direct data on the total 
charged multiplicity at the Z° pole and the corresponding number for the multiplicity of 3-jet 
events, see J JOI • Note that the multiplicity difference rises as Mh increases. We expect that 
such a large effect could help to discriminate between the Higgs signal and background events 
containing 6-quarks and in the analysis of the bb diffractive events at the Tevatron. 
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Appendix: Helicity amplitudes for gg — > gqq 

Here we outline the formalism used to calculate the gg — > gqq scattering amplitude discussed 
in Sectional We denote the colour indices of the incoming gluons by a, b, and of the outgoing 
gluon by c. The quarks colour indices are j, k. 

The gg — > gqq matrix element, which depends on the helicities, hi, and the 4-momenta, 
Pi, of the gluons and quarks, is given by the so-called dual expansion (see [SSj and references 
therein) 

M hi (Pi) jk = E( A aA fo A c ) jfc z(a, b, c) , (39) 

where the sum is over the non-cyclic permutations of a, b, c. The first factor has the same 
structure as if all the gluons were emitted from the quark line. The A, are the standard 
matrices of the fundamental representation of SU(3), and are normalised as follows 

Tr(A a A 6 ) = ^5 ab , (40) 

[A a , A 6 ] = ifabA c - (41) 

The colour-ordered subamplitudes, z(a, b,c,), are only functions of the kinematical variables of 
the process, i.e. the momenta and the helicities of the gluons. They may be written in terms 
of the products of the Dirac bispinors, that is in terms of the angular (and square) brackets 



(ab) = ( P -\pt) = J\2 PaPi \e i +>*, (42) 



[ab] = (p + a \Pb) = \j\2p a p b \e^\ (43) 

where 2p a pb = s a b is the square of the energy of the corresponding pair. If both 4-momenta 
have positive energy, the phase (j) a b is given by 

cos <p ab = sin <p ab = PJP^M, (44) 



PaPb Sab ^PtPb Sab 

with pf = p°i +pi, while the phase (f> a b can be calculated using the identity s a b = (ab)[ab\. 
Finally, the only non-zero J z = subamplitudes are 

z(a, b, c; h) = tg 3 (gc)(g ' c)(/c)2 . (45) 

S (qq) (qa) (ab) (be) (cq) 

Here g s is the QCD coupling (a s = gl/^n) and I denotes the quark (or antiquark) which 
has the same helicity as the outgoing gluon c. In particular, when A = A& = 1 while A c = \ q 
(Xq = —X q ) the numerator takes the form (qc) 3 (qc). The expression ()45j) is written for the case 
of the incoming gluons with positive helicities. If we change the sign of all helicities, then we 
have simultaneously to replace the (ij) brackets by the [ij] brackets. 

Note that in the formalism leading to (|45|) all the gluons are considered as incoming particles; 
that is, the energies of gluon c and both quarks are negative. In the case when one or two 
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momenta in the product (ab) have negative energy, the phase (f)^ is calculated with minus the 
momenta with negative energy, and then mi/ 2 is added to <j) a b where n is the number of negative 
momenta in the spinor product. 

It is clear from ()45|) that in the limit of a soft outgoing gluon c the cross section obtains an 
extra factor Eg since the numerator of the amplitude 

(qc)(qc)(Ic) 2 KE 2 g . 

This factor kills the soft gluon logarithm. The collinear logarithm corresponding to the kine- 
matics where gluon c goes in, say, the quark q direction may come only from the subamplitudes 
z(c, a, b; h) and z(c, b, a; h) where the factor (qc) in the denominator of (J45|) provides the collinear 
singularity. However this singularity is cancelled by the analogous term in the numerator. Thus 
the most dangerous background configuration (where the gluon is very hard to separate from 
the quark jet) is not enhanced by any large logarithm. 

It is straightforward to derive the analytical expression for the experimentally important 
kinematic configuration where the gluon and quark directions are aligned. This involved keeping 
just the two subamplitudes, z(c, a, b; h) and z(c, b, a; h), with the quark helicity opposite to the 
gluon c helicity (in order to have a larger numerator in (jHH)). This gives 

£ |M| 2 (J 2 = 0; colour singlet) = (46) 
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